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ABSTRACT 

We calculate a new suite of albedo models for close-in extrasolar giant planets and compare with 
the recent stringent upper limit for HD 209458b of Rowe et al. using MOST. We find that all 
models without scattering clouds are consistent with this optical limit. We explore the dependence 
on wavelength and waveband, metallicity, the degree of heat redistribution, and the possible presence 
of thermal inversions and find a rich diversity of behaviors. Measurements of transiting extrasolar 
giant planets (EGPs) at short wavelengths by MOST, Kepler, and CoRoT, as well as by proposed 
dedicated multi-band missions, can complement measurements in the near- and mid-IR using Spitzer 
and JWST. Collectively, such measurements can help determine metallicity, compositions, atmospheric 
temperatures, and the cause of thermal inversions (when they arise) for EGPs with a broad range of 
radii, masses, degrees of stellar insolation, and ages. With this paper, we reappraise and highlight the 
diagnostic potential of albedo measurements of hot EGPs shortward of ^1.3 /xm. 
Subject headings: stars: individual (HD 209458) — (stars:) planetary systems — planets and satellites: 
general 



1. INTRODUCTION 

Optical albedo measurements of solar system objects 
such as planets, asteroids, and moons have a strong tradi- 
tion in planetary science. The associated reflectances and 
their wavelength dependence have been used to deter- 
mine surface compositions and, in concert with infrared 
(IR) photometric measurements and orbital distances, 
can be used to determine radii. Orbital measurements 
can sometimes reveal masses and, given the radius, an 
object's average density can be calculated. This density 
then points to its bulk composition. In this way, optical 
albedo studies have often been central to constraining the 
interior and surface properties of solar-system bodies. 

With the discovery of extrasolar giant planets (EGPs) , 
this same program of sensing and characterization can be 
envisioned and numerous pilot programs to measure their 
optical albedos have been inaugurated (Collier-Cameron 
et al. 2002; Leigh et al. 2003; Leigh, Collier-Cameron, 
& Guillot 2003; Rowe et al. 2006). However, unless the 
EGP can be imaged separately from its primary star, one 
must contend with the severe planet— star flux ratios in 
the optical. At the distance of Jupiter from its primary, 
this ratio is ~10 -9 . Even at the distance of the closest 
transiting EGPs (^0.02 A.U.), the optical planet— star 
flux ratio is around 10 -5 (Sudarsky, Burrows, & Pinto 
2000; Marley et al. 1999). These small numbers are the 
primary reason direct measurements in the near- to mid- 
infrared using Spitzer have been more productive of late 
in the study of close-in and transiting EGPs. The cor- 
responding ratios of the absorbed optical radiation (re- 
processed by the planet into the infrared) to the stellar 
IR can be as high as 10 -4 to 10~ 2 . With such favorable 
contrasts, and using the Spitzer IRAC and MIPS IR pho- 
tometers and the IRS spectrometer, planets outside the 
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solar system have been directly measured and character- 
ized for the first time (Charbonneau et al. 2005; Deming 
etal. 2005,2006,2007; Harrington etal. 2006; Harrington 
et al. 2007; Knutson et al. 2007; Grillmair et al. 2007; 
Richardson et al. 2007; Swain et al. 2008). Researchers 
have measured atmospheric temperatures, thermal inver- 
sions, and the possible presence of water, methane and 
carbon monoxide. They have also discerned atmospheric 
features and inferred day-night climatological contrasts 
(Harrington et al. 2006; Knutson et al. 2007). 

Nevertheless, the optical offers opportunities for mea- 
surements that can complement those in the IR to pro- 
vide firmer constraints on the atmospheric properties 
of close-in EGPs. The albedo in the blue can reveal 
Rayleigh scattering off atmospheric H2 and He to deter- 
mine particle densities. Measurements in the resonance 
lines of the alkali metals Na (~0.589 /mi) and K (^0.77 
lira) directly test for their presence and are a counter 
to transit radius measurements that have already indi- 
cated the presence of Na in HD 209458b (Charbonneau et 
al. 2002). Importantly, albedo/reflection measurements 
probe different regions of the atmosphere than transit 
measurements . 

The depth of the water feature around ^0.94 /j,m ad- 
dresses both the abundance of oxygen and of water. In 
these hot atmospheres, oxygen can be mostly in CO and 
the CO/H2O ratio might someday provide a signature of 
non-equilibrium chemistry and overall metallicity. The 
latter addresses planetary origins and formation. 

The thermal inversions recently invoked to explain the 
photometric reversals and flux enhancements seen in the 
IRAC, IRS, and/or MIPS data for HD 209458b (Hubeny, 
Burrows, & Sudarsky 2003; Knutson et al. 2008; Bur- 
rows et al. 2007) and for HD 149026b (Harrington et al. 
2007; Burrows, Budaj, & Hubeny 2008; Fortney et al. 
2006) might be due to an "extra absorber" in the opti- 
cal. Such an absorber would have a strong impact on the 
optical reflection spectrum. If it were TiO, VO, and/or a 
chemical product of photolysis, such as polyacetylenes or 
"tholins," distinctive absorption features would in prin- 
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ciple be observable 3 . If the temperature inversion were 
very strong, elevating the temperature of the upper at- 
mosphere above the "plateau" value at depth, emission in 
the optical might trump absorption (Hubeny, Burrows, 
& Sudarsky 2003; Seager & Lopez-Morales 2007). On 
the other hand, if the inversion is due not to an optical 
absorber, but to a mechanical heat source like the break- 
ing of inertial or tidal waves, the resultant absence of a 
distinctive optical absorption feature in the albedo spec- 
trum would be telling. These speculations all serve to 
emphasize the fact that theory is still in its earliest for- 
mative stages, for though inversions in transiting EGPs 
were anticipated, a thermal inversion specifically in HD 
209458b, and the magnitude of that inversion, were not. 

In addition, there could be thick cloud layers of sili- 
cates and/or iron. Such clouds are seen in L dwarf at- 
mospheres and might emerge in an irradiated planet's 
atmosphere when their condensation curves intersect the 
planet's T/P profile. This could happen for the most 
severely irradiated EGPs. If these clouds form, they do 
so first near the sub-stellar point and and equatorial re- 
gions. Reflecting clouds could increase the optical albedo 
by a factor of ~2 to ~5. 

Moreover, the redistribution of heat by zonal winds 
from the dayside to the nightside lowers the dayside 
temperatures and, thereby, alters the molecular compo- 
sition. If we can measure this effect, parametrized by 
P„ (0 <P„< 0.5) (Burrows et al. 2005, 2006, 2007), we 
might have a handle on the planet's global climate. In 
the optical and very-near-IR, the variation with P n of 
the planet's flux and albedo can be a factor of 50%, de- 
pending upon wavelength. Note that the magnitude of 
the dependence on P„ of the planet— star flux ratio is 
no weaker in the optical than in the mid-IR probed by 
Spitzer. 

Therefore, due to the possible features, effects, and 
constituents mentioned above, the monochromatic op- 
tical albedos from ^0.35 to ~1.0 /Ltm for the family of 
known transiting EGPs could vary by about an order 
of magnitude. Moreover, due to physical differences be- 
tween the planets, the slopes of the albedo spectra at var- 
ious wavelengths and the contrasts in and out of various 
bands can differ by many factors. The flux at the sub- 
stellar point (hereafter "sub-stellar flux") on the known 
transiting planets ranges by ~two orders of magnitude 
from 0.03 x 10 9 ergs cm" 2 s" 1 for GJ 436b to 5.9 x 10 9 
ergs cm~ 2 s _1 for OGLE-TR-56b and the stellar types 
of the parents range from M, through K and G, to F. 
Therefore, the irradiation spectra vary substantially and 
this variety should translate into a variety of albedo spec- 
tra. In sum, optical albedo measurements complement 
near- and mid-IR measurements to more strongly con- 
strain the physical and chemical properties of tranisting 
EGP atmospheres than either can do alone. 

Among those transiting EGPs that orbit bright and/or 
nearby stars, HD 189733b, HD209458b, TrES-1, HAT-P- 
2b, GJ 436b, XO-2b, HAT-P-1, HD149026b, and TrES-4 
collectively constitute a heterogeneous class of objects 
that must have disparate albedo and flux signatures that 

3 Note that Jupiter's albedo spectrum in the blue is lowered by 
about a factor of two shortward of ~0.5 )im due to the presence of 
non-equilibrium upper-atmosphere absorbing species of unknown 
origin. 



will enable a fruitful photometric study of the atmo- 
spheres of close-in EGPs. For a non-transiting close-in 
EGP (such as 51 Peg b, v And b, or r Boo b), the ambi- 
guity in the orbital inclination and radius renders inter- 
pretation more difficult, but no less interesting. 

Recently, using the MOST (Microvariability & Oscil- 
lations of STars) micro-satellite (Walker et al. 2003; 
Matthews et al. 2004), Rowe et al. (2007) obtained a 
stringent 1-er upper bound of 8.3% on the integrated op- 
tical albedo of HD 209458b 4 . The actual measurement is 
3.8% ± 4.5%. More data on other transiting EGPs are 
sure to follow and this, along with the new HD 209458b 
constraint, make it a good time to address anew the the- 
ory of EGP albedos. In this paper, we compare the Rowe 
et al. number with a representative collection of new 
albedo models. Though we focus on HD 209458b, our 
subsidiary goal is to reignite the discussion concerning 
the influences on EGP albedos of variations in the phys- 
ical properties of their atmospheres. 

2. ALBEDO MODELS AND COMPARISION WITH THE HD 
209458B DATA 

The spherical albedo is the fraction of the stellar light 
at a given wavelength reflected into any and all angles. 
The Bond albedo is the frequency or wavelength integral 
of the spherical albedo, weighted by the incident stellar 
spectrum, and is a measure of the heat reflected. The 
Bond albedo is the quantity one uses to determine the 
planet's total energy budget. However, it is the geomet- 
ric albedo (A g ) that one measures when detecting the 
planet's flux at the Earth. Hence, this albedo is the fo- 
cus of this report and it can simply be defined using the 
formula: 

F P /F*=A 9 (^) 2 , (i) 

where F p /F* is the planet— star flux ratio at full face, 
which for transiting planets is at secondary eclipse, R p 
is the planet radius, and A is the semi-major axis of 
the planet's orbit, assumed hereafter to be circular. Of 
course, the planet— star flux ratio varies with phase an- 
gle via the phase function (Burrows, Sudarsky, & Hubeny 
2004; Sudarsky, Burrows, & Hubeny 2005), but at supe- 
rior conjunction the phase function is set to one, thereby 
fixing the definition of A g . Henceforth, we will not be 
concerned with the phase dependence and will focus on 
the full- face quantities A g and F p /F* and their wave- 
length dependence. 

Implicitly, the definition in eq. ([1]) acknowledges that 
at the high atmospheric temperatures of close-in plan- 
ets the planet itself can be self-luminous in the near- and 
mid-IR and that A g at these wavelengths can exceed one. 
This possibility is not at all anomalous, since in deter- 
mining an albedo one is taking the ratio of planet flux 
to stellar flux, at the same wavelength. In the IR, the 
planet's flux may be due more to the absorption and re- 
emission of the star's optical flux than to its IR flux and 
when the star's IR flux is the divisor (as in eq. [I}, the re- 
sulting A g can easily exceed one without violating energy 
conservation. Hence, the anomaly, if anomaly there is, is 
merely one of expectations colored by solar-system expe- 
rience, for which the bodies in question are much cooler 

4 Rowe et al. achieve a photometric accuracy better than 10 — 5 . 
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than the Sun and the concept of optical albedo is sensibly 
tied to reflection. However, for close-in EGPs, the fluxes 
are more and more thermal longward of —0.7 /im. Nev- 
ertheless, incorporating the stellar flux as a boundary 
condition and solving the equations of radiative trans- 
fer for all wavelengths self-consistently, yields values of 
Fp/F*, from which A„ as a function of wavelength can 
be derived using eq. 

To produce our models for F p /F» and A g , we em- 
ploy the self-consistent atmosphere code and solution 
techniques described in Hubeny & Lanz (1995), Bur- 
rows, Hubeny, & Sudarsky (2005), Burrows, Sudarsky, & 
Hubeny (2006), and Burrows, Budaj, & Hubeny (2008), 
to which the reader is referred for details. In Burrows, 
Budaj, & Hubeny (2008), we also describe how we now 
handle heat redistribution for a given P„ and how we 
derive the flux at secondary eclipse. In Sharp & Bur- 
rows (2007), we review our opacity libraries and chem- 
ical abundance calculations. The latter are done in the 
context of chemical equilibrium. Using these tools and 
databases, we have generated spectral models for HD 
209458b at superior conjunction as a function of metal- 
licity and P„. We have also explored the dependence on 
a possible high-altitude "extra absorber" in the optical, 
as inferred for HD 209458b by Burrows et al. (2007) us- 
ing the IRAC data of Knutson et al. (2008). Though 
the models range from 0.3 /im to 300 /im, we highlight 
in this paper only the results shortward of 1.3 /jm. 

Figure [1] portrays the albedo spectra for a variety 
of atmospheric models for HD 209458b. Superposed 
are the Rowe et al. (2007) upper limit (8.3% [l-er] = 
3.8%±4.5%), the MOST transmission curve, and the 
model bandpass-averaged albedos. Table 1 depicts these 
averaged albedos, as well as the corresponding numbers 
for a few models of HD 189733b for comparison. The 
model spectra are calculated for a given triplet of P n 
(0.1, 0.3, or 0.5), atmospheric metallicity (either solar or 
lOxsolar), and k c (0.0 or 0.1 cm 2 /g). k c is the opacity 
of a putative extra absorber which might be necessary 
to generate the thermal inversion at altitude inferred 
from the IRAC data for HD 209458b (Knutson et al. 
2008). We distribute this absorber everywhere in the at- 
mosphere at pressures lower than —0.014 bars and in the 
wavelength interval 0.43 /im < A < 1.0 /jm. This pre- 
scription should be viewed as useful, but simplistic. To 
fit the HD209458b data, Burrows et al. (2007) employed 
a constant value of K e = 0.1 cm 2 j g. 

Figure[T]incorporates our major results and manifests a 
number of features. The general trend is that the albedo 
decreases to shorter wavelengths until Rayleigh scatter- 
ing off H2 and He reverses this. This can occur near 
—0.5 /im and can be responsible for an increase in A g 
of as much as an order of magnitude from 0.55 /im to 
—0.30 /im. At its trough, A g can be as low as (or lower 
than) — 1%, but without an extra broad-band absorber 
as we have modeled it here it is only in the relatively nar- 
row alkali metal features that A g is severely low. Oth- 
erwise, the trough is around — 1% to —5%. However, 
this can correspond to values of Fp/F* less than 10~ 5 
and such numbers may be inaccessible from the ground 
for the foreseeable future. However, due to the high 
temperatures (-1500-2000 K) in HD 209458b's atmo- 
sphere (Burrows, Budaj, & Hubeny 2007), A g exceeds 
1.0 longward of —1.0 /im to — 1.1 /im. In fact, using 



the strict definition in eq. ([T]), in the mid-IR A g exceeds 
10. Note that for most of our models longward of —0.7 
/im A g is above 0.1 (10%). At —1.0 /im, it can be near 
—60%. Therefore, longward and shortward of the V and 
R bands, in the U, B, Z, Y, and z' bands, A g can be 
respectable, even high. Such high values may explain 
the very tentative 1-a lower limit to A g in the B band 
of -0.14 estimated by Berdyugina et al. (2007) for HD 
189733b (see also Table 1). ' 

However, when integrating over the MOST passband, 
our models, for the variety of values of P„ and metallicity, 
have low predicted average albedos, varying from —0.017 
to -0.07 (Table 1). These are quite "black." In fact, all 
our models are consistent with the low Rowe et al. 1-a 
upper limit of 8.3%. Reflecting clouds in the atmospheres 
of solar-system bodies, and in previous theoretical EGP 
models with clouds (e.g., Sudarsky, Burrows, & Pinto 
2000), have values of A g above 0.2 (20%), usually near 
—30-40%. Such numbers are inconsistent with the Rowe 
et al. measurement. Therefore, a straightforward con- 
clusion is that there are no reflecting clouds in the atmo- 
sphere of HD 209458b. Thin hazes are possible, but we 
estimate their optical depth in the MOST bandpass can 
be no larger than —0.1 for high single-scattering albedos 
(S = k s / (n s + n a )). This conclusion may challenge some 
otherwise reasonable solutions (Fortney et al. 2003) to 
the anomaly in the magnitude of the transit radius in the 
Na-D lines measured by Charbonneau et al. (2002). 

The Rowe et al. (2007) data provide only one number 
and the MOST transmission curve is very broad. As a 
result, we can derive little from these data about the at- 
mosphere of HD 209458b. However, the variety of models 
depicted in Fig. [T]and Table 1 allow us to make general 
predictions concerning the diagnostic potential of future 
albedo spectrum measurements from —0.3 /im to 1.3 /im 
and to discuss the dependence of A g on the physical at- 
tributes of the atmospheres of close-in EGPs. 

The relative behavior in Fig. Q] of the solar metal- 
licity models at k c = (black and two shades of gray) 
as a function of P„ demonstrates that small values of 
P„, which yield the hottest dayside temperatures, re- 
sult in higher values of A g longward of —0.65 //m. A g 
and Fp/F* increase by —30-50% as P„ shifts from 0.5 
to 0.1. This is as expected, given the thermal character 
of the longer-wavelength fluxes. However, hotter atmo- 
spheric temperatures result in higher absorption opaci- 
ties. At the shorter wavelengths, where Rayleigh scatter- 
ing is manifest, these higher absorption opacities trans- 
late into lower scattering albedos. As Sudarsky, Burrows, 
and Pinto (2000) and van de Hulst (1974) showed, A g de- 
creases quickly with decreasing single-scattering albedo, 
as much as a factor of two for a —10% decrease in S. 
We see this behavior in Fig. Q] between 0.4 and 0.5 /im, 
where the black curve (P n = 0.1) is almost a factor of two 
lower than the light gray curve (P„ = 0.5), even though 
the MOST averages are comparable. Hence, higher val- 
ues of the heat redistribution parameter result in higher 
(lower) values of A g at short (long) wavelengths. 

Comparing the brown curve (10 x solar, P ra = 0.3) with 
the dark gray curve (solar, P„ = 0.3) allows one to gauge 
the metallicity dependence. However, shortward of —1.3 
/im, the metallicity effect can be subtle. At wavelengths 
longer than —0.9 /im, higher metallicity results in slightly 
lower A g and Fp/F*, but only by —10%. Between —0.85 
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/im and ^0.55 /jm, the A g s are very similar. However, 
below ~0.5 /xm, A s s are much lower for 10 x solar metal- 
licity. The higher metallicity results in much higher ab- 
sorption cross sections and, hence, lower £. In a spectral 
region where Rayleigh scattering, not thermal emission, 
is important, as described above this translates into a 
much lower value of A g . This fact is demonstrated in 
Table 1 by comparing the B band albedos. 

Note that the transmission curves of MOST, Kepler 
(Borucki et al. 2003; Koch et al. 2004; D. Koch, pri- 
vate communication), and CoRoT penetrate progres- 
sively more and more into the red. For MOST, the 
wavelength on the red side at which transmission is 10% 
of peak is ~0.75 /im, for Kepler it is ~0.9 /im, and for 
CoRoT it is ^1.0 /im. As indicated in Table 1, the differ- 
ent predicted albedos for these different satellites may be 
crudely diagnostic of the generic wavelength dependence, 
with CoRoT sampling the longer wavelengths at which 
A g is generally higher. However, since CoRoT, MOST, 
and Kepler can't observe the same stars, until there is 
a dedicated multi-band albedo mission, we will have to 
make due with minimal constraints, however precious, 
on any particular transiting EGP albedo spectrum. 

Water absorption features reside at ^0.92-1.0 /jm and 
~1 .1-1.2 /im and are clearly seen in Fig. [TJ Their depths 
can be from ~50% to a factor of two and are clear signa- 
tures of water's presence. Nearly saturated, the metal- 
licity and abundance dependence is modest, with the 
10 x solar model lower by ^20-40% than the solar model. 
Hence, measurements in and out of these features and 
of the absolute A g s would identify water's presence and 
might constrain its abundance. 

The red (solar) and green (lOxsolar) curves in Fig. [1] 
depict two realizations of a model with an "extra ab- 
sorber" (k c = 0.1 cm 2 /g) in the optical that creates a 
thermal inversion at altitude and explains the the IRAC, 
MIPS, and IRS data for HD 209458b (Burrows et al. 
2007). While in the mid-IR these two models differ 
little (~20%), the presence of a broad optical absorp- 
tion feature can severely suppress flux where it oper- 
ates. At solar metallicity (red), the effect on A g is se- 
vere between its artificial blue (0.43 /mi) and red (1.0 
/zm) edges. However, if the metallicity is as high as 
lOxsolar, the suppressive effect of adding an extra ab- 
sorber with a given opacity is less (green curve). While 
at high metallicity the flux is generally lower in the op- 
tical than for all the models without inversions, it is a 
factor of <2 times higher than for solar-metallicity mod- 
els (red curve). This inverse dependence on metallicity 
is encouraging as a diagnostic, until one realizes that in 
an absolute sense the geometric albedos and flux ratios 
are quite low. In addition, our models place the extra 
absorber at the same pressure levels (<0.014 bars), in- 
dependent of metallicity, and, depending upon the ori- 
gin of the absorber, it is likely that its physical extent 
will somehow be metallicity-dependent. Nevertheless, 
a metallicity dependence may be a future signature of 
the metallicity/extra-absorber combination that is not so 
manifest at infrared wavelengths. Clearly, knowledge of 
the specific absorption spectrum of the extra absorber, if 
it exists, would make our predictions more specific. Such 

5 The t ransmissio n curve for CoRoT can be fo und at 
: //corotsol . obspm. f r/web-instrum/payload.param/ 



would be the case if it were TiO and/or VO (Hubeny 
et al. 2003). However, the cold-trap effect (Burrows et 
al. 2008) suggests that TiO and VO might not have the 
requisite abundance at altitude to explain the inversion 
inferred from the HD 209458b IRAC data of Knutson et 
al. (2008). This is an open question. 

There are two final features to note concerning the 
models with an "extra absorber" portrayed in Fig. [TJ 
The first is that at the shortest wavelengths the so- 
lar model has higher values of A g (red curve). These 
higher albedos are a consequence of the fact that while 
the extra absorber heats the atmosphere at low pres- 
sures (< 0.014 bars), it results in lower temperatures for 
the high-pressure plateau (Burrows, Budaj, & Hubeny 
2008, their Figure 1). This is near were the far-blue 
and near-UV photospheres reside. Lower temperatures 
result in lower absorptive opacities and, hence, higher 
scattering albedos. The result is the counter-intuitive 
enhancement in the albedos at the shortest wavelengths. 
This is less the case for the lOxsolar model, due to the 
correspondingly higher absorptive opacity. The second 
feature is that the model with an extra absorber that fits 
the HD 209458b IRAC data (Burrows et al. 2007) pre- 
dicts mid-optical albedos and corresponding planet— star 
flux ratios that are lower, not higher, than without the 
absorber, mildly at odds with the suggestion of Fortney 
et al. (2007). This is because in these extra-absorber 
models the optical photosphere is at lower temperatures 
(and pressures) than is the temperature plateau at depth, 
near where the optical photosphere would reside with- 
out the extra optical absorber. The upshot is a cooler 
emitting surface at optical wavelengths with the extra 
absorber, than without, at least for the Burrows et al. 
(2007) models that currently fit the HD 209458b IRAC 
data. Of course, a more precise spectral model for the 
thermal inversion may lead to a different conclusion, de- 
pending on the resultant T/P profile. In addition, the 
much hotter upper atmospheres expected for OGLE-TR- 
56b, TrES-4, OGLE-TR-132b, and XO-3b may indeed 
manifest this albedo enhancement (Hubeny, Burrows, & 
Sudarsky 2003; Seager & Lopez- Morales 2007). How- 
ever, our results stress caution in making qualitative 
predictions for planets with thermal inversions for the 
planet— star flux ratios in the optical. 

Moreover, if the thermal inversion is due to the break- 
ing of mechanical waves, and not an upper atmosphere 
absorber in the optical, the albedo predictions could be 
very different. In particular, A s in the mid-optical could 
be a factor of ~5-10 higher if the inversion is not due to 
a extra absorber. This difference might be a diagnostic 
of the wave breaking mechanism in planet atmospheres 
for which the near- and mid-IR data demand a thermal 
inversion. 

3. SUMMARY 

Focussing on HD 209458b, we have generated a new 
suite of albedo models for close-in extrasolar giant plan- 
ets. We find that all models without scattering clouds 
fit the Rowe et al. (2007) albedo limit for HD 209458b 
and that it is no surprise that the albedo in the MOST 
bandpass is very low. We have explored the effects of 
metallicity, of the degree of heat redistribution, and of 
enhanced optical opacities at altitude that might create 
thermal inversions such as are inferred from the IRAC 
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TABLE 1 

Band-Averaged Geometric Albedos for HD 209458b and HD 189733b 1 



Planet 




Filter 


MOST 






Kepler 




CoRoT 




B Bcsscl 






[Fc/H] 


«e\ P,i = 


0.1 0.3 


0.5 | 


0.1 


0.3 


0.5 | 


0.1 0.3 


0.5 | 


0.1 0.3 


0.5 


HD 209458b 


solar 
lOXsolar 


0.0 
0.1 
0.0 
0.1 


0.070 0.066 
0.023 
0.030 
0.012 


0.063 | 
... | 

... | 


0.111 


0.095 
0.014 
0.076 
0.022 


0.079 | 
... | 
... | 
... | 


0.136 0.117 
0.026 
0.092 
0.031 


0.097 | 
... | 
... | 
... | 


0.077 0.085 
0.057 
0.012 
0.015 


0.094 


HD 189733b 


solar 
lOXsolar 


0.0 
0.1 
0.0 
0.1 


0.044 0.044 
0.022 


0.048 | 
... | 

■ 

... | 


0.047 


0.045 
0.007 


0.044 | 
... | 
... | 
... | 


0.063 0.059 
0.017 


0.058 | 
... | 
... | 
... | 


0.084 0.088 
0.073 


0.102 



data for HD 209458b and HD 149026b. Measurements 
of hot EGPs at short wavelengths can complement mea- 
surements in the near- and mid-IR using Spitzer and 
JWST to determine metallicity, water abundances, at- 
mospheric temperatures, and the cause of thermal inver- 
sions, when they arise. Moreover, albedos at short wave- 
lengths are strongly affected by Rayleigh scattering off 
molecular hydrogen and helium and their measurement 
can be used to determine or constrain the sum of their 
number densities. Soon, Kepler and CoRoT will mea- 
sure in slightly different wavebands the optical transits, 
the light curves, and the albedos of a large number of 
known and presently unknown transiting EGPs, thereby 
providing a wealth of new data and crude spectral di- 
agnostics below ~1.0 fj,m. These data can be compared 
with models calculated using the formalism we applied 
here for HD 209458b to constrain atmospheric param- 
eters of close-in giants with a broad range of gravities, 
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Fig. 1. — The logarithm base 10 of theoretical albedos (left axis) and F p /F* ratios (right axis) as a function of wavelength (in microns) 
from 0.35 (im to 1.3 /im for a collection of models for HD 209458b. Here F s t ar = F*. Each model spectrum is for a given representative 
combination of P n (either 0.1, 0.3, or 0.5), atmospheric metallicity (either solar or lOxsolar), and "extra absorber" opacity (k c , either 0.0 
or 0.1 era? /g). The legend at the bottom right provides the matrix of models and the colors of the corresponding lines. The prescription 
for incorporating the extra absorber, when present, is described in Burrows, Budaj, & Hubeny (2008). k c is a phenomenological parameter 
with which to numerically generate a thermal inversion at low pressures. Such an inversion is inferred for HD 209458b by Burrows ct al. 
(2007) from an analysis of its IRAC data (Knutson et al. 2008) and may be due to photolytic products, TiO/VO, or some other origin. 
Superposed in blue is the MOST transmission function. Since its transmission function is so broad, ranging from ^0.35 fim to 0.8 fim, we 
provide with horizontal lines the corresponding predicted average albedos as they would be measured by MOST. The Rowe et al. (2007) 
data (3.8% ± 4.5%) are superposed in orange (dashed), with the topmost extent indicating the l-cr upper limit (8.3%). See the text for a 
discussion of the models and their implications. 



